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We study the relation between the X-ray luminosity of compact sources and the SFR of the host galaxy. Our sample 
includes 38 galaxies for which a uniform set of X-ray, infra-red and ultraviolet data from Chandra, Spitzer and GALEX 
has been collected. Our primary goals are (i) to obtain a more accurate calibration of the Lx — SFR relation and (ii) to 
understand the origin of the dispersion in the Lx — SFR relation observed in previous studies. Preliminary results of this 
project are reported below. 
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1 Introduction 

The recent star formation rate (SFR) of a galaxy is mea- 
sured from the light emitted by massive stars: given their 
short lifetimes, their luminosity is, to the first approxima- 
tion, proportional to the rate at which they are currently be- 
ing formed. The conversions between luminosities at dif- 
ferent wavelengths and SFR are usually derived from stel- 
lar po pulation synthesis models (see review by iKennicutt , 
1998). Most of the SFR indicators are known to be subject 
to various uncertainties, one of which is correction for the 
interstellar extinction. 

Another method of SFR determination is based on the 
X-ray emission of a galaxy. As X-rays are able to penetrate 
large columns of gas and dust, such an estimator is expected 
to be independent of this uncertainty of the more conven- 
tional indicators. The X-ray emission of a normal gas-poor 
galaxy is dominated by the coll ective emission o f its X-ray 
binary (XRB) population (e.g. Fabbianol 2006, and refer- 
ences therein). The latter can be divided into low- and high- 
mass X-ray binaries according to the mass of the donor star. 
The difference in mass results in different evolution time- 
scales. The characteristic time-scale of luminous high-mass 
X-ray binaries (HMXB) does not exceed ~ few x 10 7 yr, 
making them a potentially good tracer of the recent star for- 
mation activity in the host galaxy. 

Compact stellar mass X-ray sources more luminous 
than ~ 10 39 erg/s are nam ed ULXs (ultra-luminous X-ray 
sources, Mushotzkv , 2004}) . They are predom inantly located 



in star-forming galaxies. iGrimm et alj (120031) showed that 
ULXs form the bright end extension of the luminosity func- 
tion of compact sources in star-forming galaxies, the fainter 
end of which is populated by "normal" HMXBs. Due to rel- 



atively shallow slope of the HMXB XLF, the collective X- 
ray luminosity of compact sources in a galaxy is dominated 
by its brightest sources, thus ULXs drive the correlation be- 
tween X-ray luminosity and the SFR of the host galaxy. 

In the past decade several authors have discussed 
the X-ray emission of compact sour ces as an es t imator 
of the SFR of the host galaxy (e.g. 
Persic & R ephaeli. l2007t Hanalli et al 



Grimm etal., 2003 



2003). Their work 
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has established the general correlation between the X-ray 
luminosity of the galaxy (and/or of its binary population) 
and its SFR. However some discrepancy between the scale 
factors derived by different authors remained. Moreover, a 
rather large scatter in the Lx/SFR ratio has been seen in the 
individual samples. One of the limiting factors of previous 
studies was the heterogeneous nature of the X-ray and SFR 
data. With the large numbers of galaxies available nowadays 
in Chandra, Spitzer and GALEX public data archives it is 
became possible to construct large samples of galaxies for 
which homogeneous sets of multiwavelength measurements 
are available. This motivated us to revisit the problem of the 
SFR relation in order to (i) obtain its more accurate 
calibration and (ii) understand the origin of the dispersion 
seen in previous studies. 

2 The sample 

We started to build the sample of star-forming galaxies from 
the entire "Normal Galaxies" section of the Chandra Data 
Archive. It includes over 300 galaxies of all morphological 
types. We based our selection on the following four criteria. 

z) Hubble type: we selected only late-type galaxies, i.e. 
those having numerical index of stage along the Hubble se- 
quence T > 0; 

ii) Specific SFR: as star-forming galaxies contain also 
LMXBs, unrelated to the recent star-formation activity, the 
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contribution of the latter will contaminate the Lx— SFR 
relation. Although separation of the high- and low- mass 
XRBs in external galaxies is virtually impossible, the 
LMXBs contribution can be controlled in a statistical way. 
To this end we use the fact that their population is propor- 
tional to the stellar mass M* of the host galaxy (Gilfanov, 
2004 ). whereas HMXBs scale with SFR dGrimm et all 
2003), therefore the specific SFR (SFR/M*) can serve as 
a measure of their relative contributions to the X-ray lumi- 
nosity of the galaxy. Based on the scale-factors for HMXB 
and LMXB populations we estimated that a threshold of 
SFR/A/* > 1 x 10~ 10 yr _1 would be a reasonable compro- 
mise, sufficient for our purpose. In applying this criterium 
we estimated the M* and SFR from 2MASS and IRAS cat- 
alogues respectively. 

Hi) Number of compact sources: in order to perform the 
spatial analysis of the detected sources (see Sect. 3), we re- 
quired that galaxies in our sample have 10 or more compact 
sources within the isophotal diameter D25. After applying 
this criterium, all remaining galaxies had X-ray point source 
detection sensitivity of Li; m < 5 x 10 37 erg s _1 ; 

iv) Distance: for the galaxies passed through the above 
criteria, a number of redshift independent distances, deter- 
mined with different methods, were collected from NED. 



Following the results of iJacobv et alj ([1992), we chose the 
most accurate among them. For our primary sample we 
chose a distance threshold of 40 Mpc, in order to spatially 
resolve the X-ray compact sources with Chandra. This also 
allows us to identify and separate the central AGN from 
the rest of detected point sources. The thus defined primary 
sample includes 29 galaxies, with SFR values ranging from 
0.1 to 15 M Q yr- 1 in the 0.5 - 8 keV band. 

High-SFR sample. In order to explore the Lx — SFR 
relation at higher SFR, we relaxed the distance constraint 
and constructed the secondary high-SFR sample of more 
distant galaxies. A total number of 8 galaxies were found 
beyond the distance of D > 40 Mpc with the SFR in 
the 13 — 412 M & yr -1 range. They were treated as spa- 
tially unresolved in our analysis (see section 4 for details 
and caveats). Five of them are classified as ultra-luminous 
infrared galaxies (ULIRGs) and two as luminous infrared 
galaxies (LIRGs). Although they were contained in the nor- 
mal galaxies section of the Chandra archive, we verified 
that the contribution of AGN to the bolome t ric l uminosity 
in these galaxies is < 30% ( Nardini et all l2009h . to en- 
sure that the IR-based SFR measurement is not compro- 
mised by the latter. The AGN contribution in the X-ray 
band was constrained directly from the Chandra images 
(these galaxies are partly resolved by Chandra, the small- 
est one having angular size of D25 ~ 24 arcsec). Finally, 
we included 8 normal galaxies from the Chandra Deep 



Field North dBrandt et all 120011) . The K-corrected X-ray 



3 Spatial analysis 

In order to identify regions of galaxies where the HMXBs 
dominate the compact X-ray source population, we an- 
alyzed the spatial distribution of the detected compact 
sources in resolved galaxies of the primary sample. 

To minimize the contribution of CXB sources we con- 
structed cumulative and differential radial profiles of the de- 
tected X-ray point sources. These were compared with the 
surface density of CXB so u rces p redicted from log A'' 



logs' of iGeorgakakis et alj (2008). Based on this analy 



sis we chose regions where the predicted number of CXB 
sources did not exceed 30% of the total number of detected 
sources. 

Although only star-forming galaxies with SFR/M* > 
1 x 10~ 10 yr -1 have been selected, we further minimized 
the contribution of LMXBs to the total X-ray luminosity by 
locating the bulge regions of spiral galaxies and excluding 
them from X-ray and multiwavelength analysis. 

4 Multiwavelength analysis and SFR 
estimation 

The X-ray luminosity for galaxies in our primary sam- 
ple was measured as the sum of luminosities of point 
sources detected above the corresponding sensitivity limit. 
As Chandra observations of different galaxies achieved dif- 
ferent sensitivity, these luminosities need to be transformed 
to the same threshold. The luminosity of "missing" unre- 
solved sources was computed by integrating the luminos- 
ity function which d ifferential slope was assumed to be 1.6 



(Gri mm et all l2003). The normalization was determined for 



and radio lum i nosit ies for these galaxies were taken from 
Gilfanov et all (120041) . Thus, our high-SFR sample included 
16 galaxies. 



each galaxy individually, from the number of sources de- 
tected above its sensitivity limit. We chose Lx = 10 36 erg/s 
as the common luminosity threshold, i.e. the luminosities 
presented below are an estimate of the combined luminosity 
of compact sources brighter that this threshold. The count- 
to-erg conversion was done assuming a power law spectrum 
with slope of 2.0 and absorption of uh = 3x 10 21 cm -2 . 

The galaxies from the secondary high-SFR sample were 
treated as unresolved objects and their luminosity was mea- 
sured as total luminosity inside the D25 ellipse. The counts- 
to-ergs conversion in this case was determined based on the 
modeling of their X-ray spectra. 

For the purpose of studying the Lx — SFR relation for X- 
ray binaries we have subtracted from the total luminosity of 
point sources in the primary sample the contribution of the 
CXB sources. The latter wa s estimated statistically, based 
on the CXB log N - log S o flGeorgakakis et alj d2008l) . The 
resulting correction does not exceed ~ 35%. This correction 
is insignificant for unresolved galaxies from the high-SFR 
sample. 

We used data from Spitzer, GALEX and 2MASS 
archives to obtain accurate and homogeneous estimation of 
SFR and stellar mass. For resolved galaxies we measured 
fluxes of the FIR (70 ^m and 24 fim) , NIR (K s band) and 
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Fig. 1 Left: The SFR — M* plane. Different types of galaxies are plotted with different symbols, together with the lines of constant 
stellar-mass-to-SFR ratio (dashed lines). Right: The Lx — SFR relation. The solid line shows the linear relation given by the equation in 
Sec.5. 



UV (FUV and NUV) within the regions selected as de- 
scribed in Sect. 3. For unresolved galaxies from the high- 
SFR sample we measured fluxes within the D25 ellipse. 
We computed the SFR of ea ch galaxy combining 



IR- and UV- based estimates (see IHirashita et all 2003 



Iglesias-Paramo et al.Ll2006l) . This method accounts for the 
light emitted by the dust heated by young stars as well as the 
UV photons escaping the dust clouds. It has the advantage 
of being free of the model dependent attenuation correction 
of the UV emission. We also included a correction for the 
IR emission from old stellar population: 



SFR = SFR°j V + (1 - 77) SFR 



in 



(1) 



where SFR^y and SFRir are computed using standard 



calibr ations for the Salpeter IMF from 0.1 to 100 M Q 
(lKennicuttL[l998l) : 



SFR m (M yr- 1 ) = 4.6 x lO-^Lir^ergs- 1 ) (2) 
SFR° v (M yr- 1 ) = 1.2 X lO-^LNUV.ob^ergs" 1 ) (3) 

-^Nuv.obs is the luminosity at 2312 °A un-corrected for dust 
attenuation, Lm is the 8 — 1000 /im luminosity, r\ is the 
cirrus correction, i.e. the fractional contribution of the old 
stellar population to Lir. We used standard values of 77 = 
for starbursts and 7/ = 0.4 for disk galaxies. 

The SFR for HDFN galaxies were com puted based on 
their 1.4 GHz luminosities and calibration of lBelll d2003l) : 

SFR(M yr- 1 ) - 5.52 x l(T 29 Li.4GHz(erg s^ 1 ) (4) 

5 L x - SFR relation 

The obtained Lx — SFR relation is shown in the right-hand 
panel of Fig.[T] An important caveat is that we have plotted 



galaxies from both resolved and high-SFR samples in the 
same plot, ignoring the fact that X-ray luminosity of the lat- 
ter includes contribution of unresolved emission (Sect. 4), 
which potentially may be significant in star-forming galax- 
ies. This issue will be fully addressed in the future publi- 
cation. Our preliminary estimates show that contribution of 
unresolved emission does not change the X-ray luminosity 
of star-forming galaxies by a large factor. This is consis- 
tent with the fact that galaxies from the unresolved sample 
appear to lie on the extension of the relation for resolved 
galaxies. 

It has been shown that Lx— SFR relation for discrete 
sources is subject to effects of statistic of small numbers 
whic h tend to steepen the apparent slope of its low-SFR 
part ( Gilfanov et al. . 20041) . Therefore a fit to the relation 
in the entire SFR range is likely to produce biased results, 
depending on the distribution of the galaxies of the sam- 
ple across SFR. We defer detailed consideration of this ef- 
fect to the forthcoming publication and plot for the ref- 
erence in the FigQ] a linear relation Lx(ergs _1 ) = 3 x 
10 39 SFR(M Q yr- 1 ) 

6 Dispersion around the X x — SFR relation 

Despite of the special effort to create a sample with a ho- 
mogeneous set of multi-wavelength measurements and to 
minimize contamination by CXB sources and LMXBs, the 
resulting Lx— SFR relation bears a rather large dispersion 
of ~ 0.4 dex rms. In order to investigate the origin of this 
dispersion we searched for correlation of the Lx/SFR ratio 
with different quantities characterizing importance of vari- 
ous potential contaminating factors. These included the sen- 
sitivity limit of Chandra data L\\ m , contribution of CXB 
and LMXB sources characterized by the ratio of the stud- 
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Fig. 2 Search for correlation of Lx /SFR ratio with various potential contaminating factors: variations in sensitivity limit 
of X-ray data (a); contribution of CXB sources (b) and LMXBs (c); dust attenuation (d) and inclination of the galaxy (e). 
Resolved galaxies are shown by squares, unresolved - with triangles. See Sect. 6 for details. 



ied area and enclosed stellar mass to the SFR, dust attenua- 
tion characterized by the Ltr/£tjv ra ti° an d inclination of 
the galaxy i. The results are shown in Fig. [2] In these plots, 
the symbols without error bars represent individual galax- 
ies, and symbols with error bars show binned values. The 
vertical error bars represent the bin rms. We searched for 
statistically significant correlations in these data using the 
Spearman's rank correlation test and found none, with the 
possible exception of the dust attenuation where the associ- 
ated probability was on the ~ 0.4 — 6% level, depending on 
the samples combination. As it is obvious from the plot, the 
correlation is such that galaxies with stronger dust attenua- 
tion have smaller Lx/SFR ratio. This is to be expected. We 
note however, that in order to suppress the 0.5-8 keV lumi- 
nosity by a factor of 10, a column density of ~ 5 x 10 23 
cm -2 is needed. 

7 Conclusions 

We presented preliminary results of the study of the rela- 
tion between SFR and X-ray luminosity of galaxies based 
on Chandra, Spitzer, GALEX and 2MASS observations of 
nearby star-forming galaxies. We found an Lx — SFR re- 
lation that is in a general agreement with results of earlier 
studies. Despite our attempt to minimize the contribution 
of various contaminating factors, we still see the relatively 
large dispersion in the Lx — SFR relation, of the order of 
~ 0.4 dex rms. We investigated the possible role of var- 
ious effects and concluded that the observed dispersion is 
unlikely to be caused by any of the obvious contaminating 



factors and may have a physical origin. Plausible candidates 
are differences in the recent star-formation history, metallic- 
ity, attenuation by the dust and gas etc. 
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